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Abstract. - We determine analytically how Efimov trimer states are modified by three-body losses 
within the model of Braaten and Hammer. We find a regime where the energies approach the 
positive real axis and the decay rates vanish. 



Introduction. — Motivated by nuclear physics, Efimov discovered that three particles 
with short-ranged two-body interactions of scattering length \a\ = 00 can support an infinite 
series of weakly bound trimer states [1]. In contrast to systems of nucleons or of 4 He atoms, 
ultracold gases of alkali metal atoms offer the possibility to tune a by using a Feshbach 
resonance, and thus to go more deeply into the limit \a\ — > 00 where the Efimov effect sets 
in [2-9]. However the situation for alkali metal atoms is complicated by the existence of 
deeply bound dimer states: a weakly bound trimer necessarily decays into a deeply bound 
dimer and a free atom [6, 10]. This three-body loss process has prevented so far to produce 
weakly bound trimers experimentally The first evidence for the existence of a weakly bound 
trimer state is the recent observation of a peak in the three-body loss rate from an ultracold 
atomic Bose gas at large negative a [2]. Such peaks were predicted to occur for the values 
of a where the energy of an Efimov trimer reaches zero [5,11]. More recently, a similar 
effect was studied experimentally and theoretically in a three-component Fermi gas [12-15]. 
The simplest description of three-body losses is the model of Braaten and Hammer, where 
three incoming low-energy unbound atoms are reflected back as three unbound atoms with 
a probability e~ 4?? * and recombine to a deeply bound dimer and an atom with a probability 
1 — e n * , rj* being the so-called inelasticity parameter, which depends on the short-range 
details of the interaction potential and of its deeply bound states [3,4, 11]. This model was 
used to obtain the decay rate of Efimov trimers in the regime of small inelasticity parameter 
77* <C 1 [4, 10], as well as 3- and 4-body scattering properties [3,4, 16] and the decay rate 
of efimovian 3-body states in a harmonic trap [17,18]. This model is expected to become 
exact in the zero-range limit where |a| and the inverse relative momenta between atoms 
in the initial state are much larger than the range and effective range of the interaction 
potential [3,4,10,11,19]. 

In this Letter, we determine how Efimov trimers are modified for an arbitrary inelasticity 
parameter 77* by solving analytically the model of Braaten and Hammer. We find that the 
Efimov spectrum rotates counterclockwise in the complex plane by an angle proportional 
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to 77*. When ry* reaches the critical value where this angle equals it, the discrete states 
disappear into the continuum. When rj* approaches this critical value from below, the 
energies approach the positive real axis, so that the decay rates tend to zero. Thus a large 
inelasticity parameter 77* ~ 1 can paradoxically give rise to long-lived three-body states. 




Fig. 1: Without three-body losses the inelasticity parameter 77* = and the energy of an Efimov 
state lies on the negative real axis. The effect of three-body losses is to rotate this energy in the 
complex plane by an angle 2r?*/|s|, where \s\ = 1. 00624.... When this angle approaches n, the 
decay rate T = — 2lmE/h vanishes. 



Efimov states without three-body losses. — We start by reviewing Efimov's solu- 
tion of the three-body problem in the absence of three-body losses. We restrict for simplicity 
to the case of three identical bosonic particles, where the wavefunction i/>(ri, r2, r^) is com- 
pletely symmetric. 1 In the limit where the range of the interaction potential is negligible 
compared to \a\ and to the typical de Broglie wavelength, the interaction can be replaced ex- 
actly by the zero-range model, see e.g. [1,3,7,9,16,17,21-24]. An cigenstate of the zero-range 
model solves (i) the Schrodingcr equation 

i— 1 

when all interparticlc distances r,j are different from zero, and (ii) the Bethe-Peicrls bound- 
ary condition, imposing that there exists a function A such that 

^(rx.ra.ra) = ( — J A(R y -,r fe ) + 0(r y ) (2) 

in the limit — where particles i and j approach each other while the position of their 
center of mass and of the third particle are fixed. In what follows we assume that 

The case of different masses and statistics can be included without difficulty: only the value of s is 
modified in the hyperradial problem [20]. 
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\a\ is much larger than the typical distance between particles, so that we can set \a\ = oo. 
Equations (1,2) are then solved by Efimov's Ansatz [1] 



^(ri,r 2! r 3 ) =F(R) [1 + P 



13 







A3) — sili 


s arctan (J^j 


/ rp 





(3) 



where Py exchanges particles i and j, s ~ i ■ 1.00624 is the only solution s € i ■ (0; +00) of 
the equation scos(s7r/2) — 8/\/3 sin(,s7r/6) = 0, the Jacobi coordinates are r = ||r2 — ri|| 
and p = |[2r*3 — ri — 1*2 1| / a/3, the hyperradius is R = \J (r 2 + p 2 )/2, and the hyperradial 
wavefunction F(R) solves the hyperradial Schrodingcr equation 



d 2 
dR 2 



R dR 



s 



F(R) = ^EF(R). 



(4) 



In the limit R — > where all three particles approach each other, the attractive effective 
potential s 2 / R 2 diverges strongly, and it is necessary to impose a boundary condition on the 
hyperradial wavefunction F(R), as first realized by Danilov [25]. This boundary condition 
is conveniently expressed as 



F(R) oc ( — 
V ; R^o \R t 



(5) 



where the three-body parameter R t depends on short-range physics and is a parameter of 
the zero-range model. The solution of Eqs. (4,5) is given by the famous geometric series of 
Efimov trimers 



2h 2 

mR 2 



exp 



r ar g r(l + ,s) 



exp n 



2- 



with a (here unnormalizcd) wavefunction 

F{R)=K s (kR) 
where K is a Besscl function and n is defined by 

h 2 k 2 



E 



2m 



(6) 



(7) 



(8) 



with the determination k > ensuring that F(R) decays exponentially for R — > 00. Efi- 
mov's spectrum (6) is unbounded from below, in agreement with the Thomas effect [26] 
and with the fact that the spectrum for an interaction of finite range b coincides with 
Efimov's spectrum only in the zero-range limit, i.e. for weakly bound trimers satisfying 
\E\ < h 2 /(mb 2 ) [1,3,4,21,27-29]. 

Effect of three-body losses. — We now determine how Efimov's result is modified by 
three-body losses within the model of Braaten and Hammer. The only difference between 
the Braatcn-Hammer model and the zero-range model is that the boundary condition (5) is 
replaced by [3,4] 



R 



Rt 



where the inelasticity parameter 77* and the three-body parameter R t are parameters of the 
Braatcn-Hammer model, whose values depend on the details of the finite-range interactions 
which one wishes to model. 2 The physical meaning of (9) is that the outgoing wave cx R s 



2 E.g. for 133 Cs atoms near the —11 G Feshbach resonance, the fit of the theoretical result of [11] to the 
experimental data of [2] gives r/„ ~ 0.06 [2] and Rt ~ 30 nm [2, 17]. 
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has an amplitude which is smaller than the amplitude of the ingoing wave oc R~ s by a 
factor e -2 ''* , i.e. three ingoing atoms are reflected with a probability e~ 4 ''* and are lost by 
three-body recombination with a probability 1 — e~ 4?) * . 

The Braaten-Hammer model is expected to become exact in the zero-range limit [3, 
4, 10, 11]. This is supported by numerical calculations of the three-body loss rate from 
an atomic gas with finite-range interaction potentials, which are in good agreement with 
the prediction of the Braaten-Hammer model in the zero-range regime [19]. Moreover this 
can be explained using the adiabatic hyperspherical description [3,6,30,31]: in addition 
to the 'atomic' adiabatic channel responsible for the Efimov effect, there is one 'molecular' 
channel associated with each deep two-body bound state; one can thus look for decaying 
Gamov states with a complex energy E by imposing outgoing boundary conditions in the 
molecular channels; the coupling between the atomic channel and the molecular channels 
is only effective at distances on the order of the potential range, where the wavefunction is 
insensitive to E in the zero-range limit; and matching this short-distance wavefunction to 
the atomic-channel wavefunction at hyperradii much larger than the range and much smaller 
than the typical atomic-channel de Broglie wavelength yields the boundary condition (9). 

In the absence of losses (77* = 0), the boundary condition (9) reduces to Efimov's bound- 
ary condition (5), and one recovers Efimov's spectrum (6): 



In presence of losses (77, > 0), we solve the hyperradial problem (4,9) with the additional 
boundary condition that F(R) must decay quickly enough at infinity. 3 Using the known 
properties of Bessel functions [32] , we obtain the energies 



i. e. the spectrum is rotated in the complex plane counterclockwise around the origin by an 
angle 2r7*/|s|. The result (11) only holds if the angle 2t7*/|s| < n, i. e. for 77* < ?7* c with 



E n (v* = 0) = E\ 



(10) 




(11) 



V*c = = 1.5806..., 



(12) 



while for 77* > ?/* c there is no normalisablc solution. 

The wavefunction is still given by Eqs. (7,8), now with the determination 



Rere > 



(13) 



of the sign of re, which ensures that the wavefunction is normalisablc since 



F(R) cx 



e 



(14) 



R^oo 



The decay rate 



r = -llmE 



(15) 



is given by 




(16) 




3 More precisely the normalisation condition is dRR\F(R)\ 2 < 00, since the quantity 

J dr± dr2 | VC r l > r 2> r 3)| 2 i which does not depend on r3, has to be finite. 
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The proportionality between T and the energy was first observed in numerical calculations 
with finite-range interactions and explained using the adiabatic hyperspherical description 
in [6]. _ 

An interesting effect occurs when r)* approaches the critical value r\* c from below: the 
energies approach the opposite of the energies of the Efimov states without losses 

EM — > |£fo.=0)|, (18) 

i7«-n?,c 

so that the decay rates tend to zero. Moreover, the sizes of the states diverge: since the 
imaginary part of k tends to a positive value and its real part tends to + , the behavior (14) 
of the wavefunction at large R is an ingoing wave with a slowly decaying envelope. Physically, 
it is not surprising that the states become increasingly delocalised before disappearing into 
the continuum. Since the wavefunction is normalized, this divergence of the size implies that 
the probability for the three particles to be close to each other vanishes, so that T vanishes. 

This effect occurs within the effective low-energy model of Braaten and Hammer. We 
thus expect that it also occurs for finite-range interactions supporting deeply bound dimers, 
provided the two-body interaction potential is tuned in such a way that 77* is slightly below 
?7* c ; this could be checked numerically using the methods of [5,6,8,19,31]. 

Experimentally, it is rather unlikely to find a Feshbach resonance close to which r?* is 
slightly below r)* c . However this regime may become accessible if interatomic interactions 
are tuned using an additional control parameter, e. g. an electric field [33]. 

* * * 
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